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DETERMINATION OF SURFACE HEAT FLUX USING A 


SINGLE IMBEDDED THERMOCOUPLE 

« 

By S. D. Williams and D. M. Curry 
Lyndon B. Johnson Space Center 


SUMMARY 


A numerical method by which data from a single embedded thermocouple can 
be used to predict the transient thermal environment for both high- and low- 
conductivity materials is described. The results of an investigation performed 
to verify the method clearly demonstrate that accurate transient surface heat- 
ing conditions can be obtained from a thermocouple 0.762 centimeter from the 
surface in a low- conductivity material. Space Shuttle Orbiter thermal protec- 
tion system materials having temperatvire- and pressure-dependent properties and 
typical Orbiter entry heating conditions were used to verify the accuracy of 
the analytical procedure . 


INTRODUCTION 


The design and development of a reusable thermal protection system (TPS) 
for the Space Shuttle is dependent on a detailed knowledge of the aerothermo- 
dynamic environment to which the TPS will be exposed. The TPS thermal perform- 
ance is normally obtained from exhaustive plasma arc and radiant heating tests 
to establish reuse temperature and thermal response characteristics . In a pre- 
vious study by C\irry and Williams (ref. l), a nonlinc-ar least squares method 
was developed for the estimation of thermal property values from experimental 
in-depth temperature data. The current investigation was motivated by the re- 
sults of the previous analysis and by the opinion that calculations of surface 
heating rates and temperatures would be of considerable value in the test 
(fli^t and ground) and development phases of the Space Shuttle. 

The calculation of surface heat flux and surface temperature from an in- 
depth temperature history measurement is called the inverse heat conduction 
problem and has been discussed by nxmerous investigators (refs. 2 to 12). An 
excellent discussion of previous investigations (refs. U to 8) for solving the 
inverse problem can also be found in reference 3. In particular, Beck and Wolf 
(ref. 2) presented a method of solution using least squares and future tempera- 
tures. In a later publication, Beok (ref. 3) presented a technique using 
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nonlinear estimation in the solution of the inverse problem. Howard (ref. 9) 
developed a numerical method for determining the heat flux to a thermally thick 
wall with variable thermal properties using a single in-depth thermocouple . 

His best results were obtained for temperatxire measurements close to the 
heated surface in conjunction with a small computing interval. 

Comette (refs. 10 and 11), in analyzing the Project Fire calorimeter 
data, developed a transient inverse solution vhat required curve fitting of 
the basic temperature data. Comette 's solution accounted for variable ma- 
terial properties and yielded a closed-form euialytical expression for the local 
surface heat flux at a given instant of time. The temperature-time data for 
several thermocouples embedded in a calorimeter plug we ’e smoothed and the data 
replaced with a polynomial equation for temperature (at a particular thermo- 
couple location) as a function of time. Imber and Khan (ref. 12) developed a 
closed-form inverse solution for consteuit properties and heat flux using two in- 
depth thermocouple readings. The solution was obtained by means of Laplace 
transform techniques in which the input thermocouple data were approximated by 
a temporal power series and a second series of error functions . 

The purpose of this report is to disseminate the results of an investiga- 
tion for predicting the surface conditions (heat rate /temperature) for a coated 
insulative material (low conductivity) and for a high-temperat\ire coated carbon 
(hi^ conductivity) material. The method, -ring a single embedded thermocouple, 
accounts for variable thermal properties (as functions of temperature and pres- 
sure) as well as for the effect of radiation losses and in-depth conduction. 

In addition, the results can be obtained with approximately the same computa- 
tional time required to solve the thermal model using a known heat rate. 

As an aid to the reader, where necessary the original units of measure 
have been converted to the equivalent value in the Systeme International 
d'Unit^s (si). The SI units are written first, and the original units are 
written parenthetically thereafter. 


SYMBOLS 


A, B, C quadratic coefficients 

a, b, c, d coefficients of square temperature matrix 
Cp specific heat of material at constant pressure 

c constant defined by equation (ll) 

f function defined by equation (lO) 

i individual measurements 

k thermal conductivity of material 
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thickness of material 
material designator 
node identifier 
point 
heat flux 

convective heating rate 
calculated convective heating rate 

actual convective heating rate 

heat rate 

temperature 

temperature at end of time step 

calculated value of temperature at node r 

'■.hermocouple 

time 

distance 

Nevton-Raphson expression 
computing time, irterval 
node thickness 
convergence tolerance 
emissivity 
arbitrary value 
density of material 
Stefan-Boltzmann constant 
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Subscripts : 


i location 

J thermocouple location 

m node Identifier 

8 surface 


THEORETICAL FORMULATION 


Background 


The heat 
standard heat 


conduction vlthln a one-dlmenslonal thermal model Is based on the 
conduction equation found in numerous heat-transfer textbooks. 
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( 1 ) 


where p is the density of the material, is the specific heat at constant 

pressure, k is the thermal conductivity, and T is the temperature of the 
material at location x in the materisd at time t. 


The solution can readily be obtained if the boundary conditions and the 
initial temperature profile are known. If a heat rate is imposed on 

the surface of a material of thickness L and if the backwall boundary is in- 
sulated, the boundary conditions are 



x»0 


”%et 


(2a) 



= 0 

x*L 


(2b) 


with the initial condition given by 


T(x,t) 


t»0 ' 


(2c) 
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If thi* conducting medium is homogeneous, it can easily be shown that the 
solution obtained from equations (l) and (2) is unique. Under these condi- 
tions, the equations are linear and analytlced techniques can be used to solve 
for the temperatures. When the thermophysical properties are temperature de- 
pendent, however, the equations become nonlinear and recourse is usually made 
to numerlced methods. An implicit n\imerlcal solution has been used in this 
investigation. 


It is assiuned that a set of thermocouples 


TCj has been placed at known 


locations x in the material. This assumption implies that a set of interior 

J 

temperature-time histories T^(Xj,t) exists. If the temperature date, are 

available at x., then by solving the boundary condition 
J 





= -q 


(3) 


for 4 together with the unknown temperatures, the heat rate at location x. 

can be found lirectly (ref. 13). If the temperature-tirae history data are 
available at the surface, then the convective heating rate can be found from 


^et 


conv 



(M 


where 

rate , 
Stefan 


q^^^ is the total heat rate input, is the convective heating 

c is the total hemispherical emittance of the surface, a is the 
-Boltzmann constant , and is the temperature at the surface . 


Because temperature data generally are not available at the surface, an 
alternate technique must be used to determine the unknown heat rate by means 
of thermocouples located within the material. A technique developed by Iraber 
and Khan (ref. 12) was successfully used as a first approximation in the de- 
termination of a constant heating rate; however, the method was unsuccessful 
for a time-dependent heating rate. 


The problem in solving for the heating rate at the surface when the tem- 
perature at the surface is unknown is that one of the required boundary con- 
ditions is unavailable. There is no difficulty in solving for temperatures at 
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If the thermocouple is located on the surface. 


the 


%et 


calculated directly from the tridiagonal matrix (appendix A). 


value can be 
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any location between two thenaocouplea because both boundary conditions are 
known. The temperature at x ■ Is given by 



(5) 


and equation (2b) can be used at x ■ L. 


New Technique 

A technique has been developed to solve for the heat rate at the surface 
at each time step rather than to solve for the entire history. The technique 
is iterative; initially, a Newton~Raphson technique is used, then a quadratic 
fit. To derive the technique, it is necessary first to examine the finite dif- 
ference approximation at the surface. 



T' - T’ , 
m m»l 


t^m t^p 

m 

2 At 




( 6 ) 


where At is the interval of conq)uting time, ' denotes the temperature value 
at time t + At, m is the node identifier, and Ax is the node thickness. 
This form of the equation is for a composite material where I is the material 
designator. 

This expression can be rearranged into the tridieigonal form 


a T' , •►bT'-t'cT'+d = 

m nr»-l mm m m-1 m 

where the coefficients of the ten:; -ecure matrix are 
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(8b) 
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TiT * Tit 
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(8c) 


. £^p 

, 1 fjn ^ , 

m " 2 At *^et 


(8d) 


The Nevton-Raphson expression is 




' 7^ 


(9) 


where the function 


f(z) 


a T' , b T' ♦ 
m m+1 m m 


c T' , 
m m-1 


+ d = 0 
m 


(10) 


and the constemt c is limited to 


0 < c < 1 


( 11 ) 


to prevent excessive corrections. 

Assuming equation (I 4 ) is valid, then equation (lO) may be rearranged as 


f = 



l^m £^p 


Ax, 


T 

m 


•f 


b 

m 




conv 



+ T' 
m 


2b At 
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(12) 



Obtaining f ■ df/d4 yields 


df 


3T’ 
3f m 


3T* 34 
n ^ 


et 


(13) 


Then, from equation (12) 


3f 

3T' 

m 


1 



arid 


3T' 

m 




1 
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Equations (12), (13), and (lU) are combined to yield 




‘•('^et) 


(luo) 


(l»4b) 


(15) 


Because the thermocouples generally are located internally rather than at 
the surface, application of the Nevton-Raphson technique will result in a mon- 
otonic approach to the desired solution, but convergence is very slow and may 
require an excessive n-jraber of iterations. Therefore, as a practical solution, 
it is necessary to switch to an alternate technique, such as a quadratic fit. 
Use of the Newton-Raphson technique will usually guarantee that the newly cal- 
culated values of q_ ^ will differ from the previous iteration; however, if 

for some reason this difference does not occur, it can be forc?d by taking 
4 q i n for the two additional val'jes required, where n is some arbitrary 

value such as 0.1. The initial value of 4^^^ usually chosen to be the 

converged value for the previous time step. For the first time step, an arbi- 
trary value such as 1 may be used. The iterative process is halted when 


T» - T»(t) < 6 T»(t) (16) 

r n - n 
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where is t) -* calculated value of the temperatixre at node r corresponding 

to the thermocouple location, temperature given by the thermo- 

couple at that time, and 6 is the relative convergence tolerance. It should 
be noted that both the net heating rate and the surface temperature have unique 
solutions, whereas the convective heating rate is dependent on the assvuned 
value for emissivity. 

NumericeLL difficulties arise in determining the surface heat rate or the 
surface temperature from data based on interior thermocouples. This difficixlty 
is partly due to the timelag imposed on the system resulting from the finite 
distance between the surface and the thermocouple location. Another factor is 
the damping of surface changes at the thermocouple location. Other errors that 
may arise are due primarily to the method used for approximating the thermal 
model, the magnitude of magnitude of the thermocouple tempera- 

tures. Of course, it is assvuned that time steps compatible with the physical 
system would be used for the thermal model. 


As previously discussed, it was necessary to switch from the Newton-Raphson 
method to a quadratic fit to insvire a converged solution. This switch is ac- 
complished by assuming three points: P^(x^,y^), t and P^(x^,y^), 


where x. = [T' - T*(t)]. and y. = 
1 r n 1 1 

formed to include all three points . 




A quadratic equation can be 


The equation is 


y 



+ Bx + C 


(IT) 


where A, B, and C are quadratic coefficients. By substituting each of the 
three points into equation (lT)i one has a system with three equations and 
three unknowns (the coefficients A, B, and C). 

After the coefficients have been determined, a point on the quadratic 
curve may be found. The points have been formed such that solution is at 
y = C or at the point Pj^(0,C). 

If, in evaluating the original function with C to obtain a new Xj^, the 

solution is not within the desired tolerance as required by equation (l6) (i.e., 
|xj^ I < e), Pj^ is substituted for one of the previous points (e.g., one with 

the largest |x|), and the coefficients are determined again. The process is 
repeated until the desired solution is obtained. Normally, this process re- 
quires four iterations since the Newton-Raphson technique was converging to 
the desired solution. 
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NUMERICAL EXPERIENCE 


The numerical method discussed in the previous section has been evaluated 
for typical Space Shuttle Orbiter materials and environments. In general, the 
Space Shuttle Orbiter reusable TPS consists of reusable carbon-carbon (RCC) 
where surface temperature exceeds 1533 K and reusable surface insulation (RSI) 
for areas with maximum svirface temperatures of less than 1533 K. The thermo- 
physiceil properties used in the investigation are presented in appendix B. 

Since the objective was to develop a program to be used in both preflight and 
postflight analyses of Space Shuttle missions, cases typical of both ground and 
flight-test environments were considered. 


The RSI Thermal Model 

The RSI thermal model (fig. l) consists of 5 .08-centimeter-thick primary 
insulation with embedded thermocouples . The insulation used for this model is 
an all silica, rigid composite material (LI-900) having a density of ll|U.l8 

3 

kg/m . The boundary conditions are assumed to be a heat rate on the surface 
and to be adiabatic on the backwall. An initial temperature of 29^ K and an 
emissi^rity of 0.8 on the surface with radiation to space were assumed for this 
me del . 

Constant-heating rate .- Groiond tests are generally characterized by both 
constant-heating and constant-pressure conditions. Therefore, this investiga- 
tion was conducted to assess the effects of the single thermocouple location. 
Four different cases were \:ised for comparison: (l) the lead thermocouple was 

on the surface, (2) the lead thermocouple was 0.25^ centimeter from the surface, 
(3) the lead thermocouple was 0,508 centimeter from the surface, and {k) the 
lead thermocouple was 0.7^2 centimeter from the surface. In each case, the 
data were generated for the lead thermocouple by losing a surface heat rate of 

136 188 W/m^. 

Approximately five-place accuracy was maintained at each time step for the 

-5 -5 -6 

first three cases using convergence criteria of 10 ,10 , and 10 , respec- 

tively. For case L, only four-place accuracy was maintained with a convergence 

criterion of 10 This loss of accuracy with a reduction in 6 clearly dem- 
onstrates the damping effect as the thermocouple is placed farther from the 
heated surface. 

Variable-heating rate .- Flight conditions combine the effects of both a 
transient heat rate and a changing pressure environment. Therefore, these 
effects are discussed vising the numerical method. .The material used for this 
investigation was LI-900 with and without a svirface coating. When the coating 
was used, it was assumed to have a uniform depth of 0.0381 centimeter. Three 
trajectories (fig. 2) were applied to the RSI: a short-range trajectory Ides- 

ignated as trajectory A), a medium-range trajectory (designated as trajectory 
b) , and a long-range trajectory (designated as trajectory C). 
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The primary concern in this investigation was to determine the effects of 
thermocouple depth x. and of the convergence criterion 6 on the accuracy of 

-h _5 

surface heating rate calculation. The values used for 6 were 10 , 10 , 

and 10 Thermocouple depths of 0.25^, 0.508, and 0 .t 62 centimeter from the 
heated surface were used for uncoated materials; for the coated material, the 
thermocouple was placed 0.0381 centimeter (just behind the coating) and 0.25^ 
centimeter from the heated surface. 

The effects of the convergence criterion and thermocouple depth on the 
2 

average error in the surface heating rate for trajectory A can be seen in 

2 

table I. The analysis was performed '^•>r a consteint pressure of 101 325 N/m 
(l atmosphere). The effects of the average error for a constant 6 are shown 
in figure 3. Basically, for each additional 0.25^ centimeter in depth of the 
thermocouple, the convergence criterion must be decreased by a factor of 10 to 

2 

maintain the same relative accuracy of approximately 90 W/m . 

The heating rate provided by trajectory B was used to show the effects of 
pressure (used in conductivity calculation) on the accuracy of the surface 
heating rate calculation. The pressure curve used with trajectory B is shown 
in figure U. The summary of the results comparing the accuracy of constajit- 

2 

pressure (lOl 325 N/m (l atmosphere)) conductivity with variable-pressure con- 
ductivity can be seen in table II. The accuracy for this trajectory was ap- 
proximately the same as was obtained for trajectory A. The differences between 
press\ire-dependent and constant-pressure results are considered small enough to 
have almost no effect. In comparing the average error figures, one must also 
conclude that the differences are of almost no consequence. 

The effect shown by coating the surface of the material demonstrates the 
feasibility of the method for calculating surface heating conditions using 
actual Space Shuttle materials . The results also indicate the type of accuracy 
anticipated if the thermocouple is placed very near the surface. 

The heating rate for trajectory C is representative of a Space Shuttle 
Orbiter location characterized by a transition to turbulent flow. No pressure 
curve was used for this trajectory since it was intended primarily to demon- 
strate the capability of the program to follow rapid changes in the heating 
rate associated with a transition from laminar to turbulent flow heating 
(fig. 2). The embedded thermocouple was located 0.25^ centimeter from the 

2 

heated surface. The average error of 1.3 W/m is an order of magnitude smaller 
than the average error for trajectories A and B, but this reduction is due to 


3 / 1 n ^ ^ ~ 

The average error equals J ~ 2^ (q. - q*) > where q. 

1 ^ i=l ^ ^ ^ 


is the calcu- 


lated convective heating rate, q* is the actual convective heating rate, and 
n is the total number of individual meas\xrements i taken. 
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the order of magnitude difference in heating rates. The data presented in 
table III reveal that the percent error is approximately the same as that 
achieved for the other trajectories . 


. The RCC Thermal Model 

The RCC thermal model (fig. 5) consists of a 0.635-centimeter-thick layer 
of RCC, an airgap, and a 5 *08-centimeter-thick layer of RSI insulation. Al- 
though this model does not properly account for the internal radiation of the 
actual Space Shuttle configuration, the thermal response is representative 6ind 
provides for a good analyticed experiment using a high-conductivity material. 
The trajectory (designated as trajectory D in fig. 6) is typical for the fuse- 
lage nose area. An emissivity of 0.85 was vised for all surfaces for radiation 
transport properties. 

Although the airgap behind the carbon was vented, no pressure curve was 
used for this trajectory because it was used primarily to demonstrate the capa- 
bility of the program to follow rapid changes in the heating rate through a 
hi^-conductivity material. The thermocouple was installed on the backwall of 
the RCC (J .635 centimeter from the heated surface). The average error of 30 
2 

W/m is of the same order of magnitude as that for trajectories A and B. The 
percent error is approximately the same as that achieved in the RSI trajector- 
ies (table IV). 

For the RSI and RCC analytical verification studies already discussed, the 
same mathematical formulation controls the thermal model for both the genera- 

3 

tion of data and the prediction of the heat flux . However, under actual test 
conditions, some of the assumptions used in formulating the mathematical model 
are violated. Examples of such variations are measurement errors associated 
with vincertainties in the location of thermocouples, void areas surroxinding 
the thermocouples, and thermophysical property deviations. The effect of any 
of these errors on the predicted heat flux was not investigated because de- 
tailed error analyses of this type can be fovind in references 9 and lU. 


CONCLUDING REMARKS 


An inverse solution technique using a single embedded thermocouple has 
been developed for predicting the transient thermal environment to which the 
Space Shuttle Orbiter thermal protection system is exposed during entry. The 
accuracy of the numerical method has been demonstrated for both a high- and a 
low-conductivity material by comparison with analytically generated data. The 


■^Comparisons with known values of heating rates indicated a maximum error 

-2 -3 

of 7.8 X 10 percent, but most values were within 1 x 10 percent. 
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procedure developed is quite general and heis been Incorporated into a pre- 
viously developed program used to compute thermal conductivity values from ex- 
perimental data. Thus, a capability now exists for computing surface condi- / 

tions (heat flux and/or temperat\ire ) and thermed. conductivity values using the 
data from a single experiment. 

In addition to being sensitive to the peak surface-temperature level, the 
reusable svirface insulation material is sensitive to the rate of change of the 
surface temperature, which induces in-depth thermal gradients and stresses 
within the reusable surface insulation tiles . The inverse program can be used 
to define thermal shock test conditions by sin^jly specifying a surface- 
temperature gradient and then solving for this required treuisient heat flux. 


Lyndon B. Johnson Space Center 

Nationeil Aeronautics and Space Administration 
Houston, Texeis , February l6, 1976 

986-15- 31-OU-T2 
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TABLE II.- COMPARISON OF AVERAGE ERRORS IN DETERMINING 


HEAT RATES FOR TRAJECTORY B 


Materied 


Lead thermocouple Convergence 
depth, cm criterion 6 


Variable pressure 


Average 
error, W/m^ 


Uncoated 

LI-900 

Coated 

LI-900 

0.25^ 

.508 

.762 

.0301 

.2^h 

10-5 

10"^ 

10-« 

lo"^ 

10 "^ 

21.6 

150 

lOU 
20. U 
31.1 


Constant 

pressure 


Uncoated 

O. 25 U 

10 "^ 

17.6 

LI-900 


c 



00 

0 

LT\ 

10“^ 

1 U 9 


.762 

10-' 

55 

Co ■'ted 

.0381 

10 "^ 

20.2 

LI -900 


c 



. 251 * 

10"^ 

28.9 


L5 


















TABLE III.- COMPARISON OF CALCULATED AND KNOWN 
HEAT RATES FOR LI -900 FOR TRAJECTORY C^ 


Error, 

2 


Known heat 
rate, w/m^ 

Calculated heat 
2 

rate, W/m 

1 759 

1 759.15 

h 857 

k 857. 

16 569 

16 569.8 

50 8U3 

50 8U2.9 

65 257 

65 252 

71 039 

71 839. 

7U 790 

7U 789.8 

73 31U 

73 313.9 

65 711 

65 709.7 

53 908 

53 906.7 

liO 516 

uo 515. 

26 557 

26 555.9 

16 116 

16 115.6 

25 762 

25 762.5 

31 550 

31 550.5 

22 130 

22 129.3 

5 777 

5 776.7 







Error, 

percent 


- 0.00306 

-.00061 

-.00170 

.OOllU 

.00721 

-.00255 

.00019 

.00075 

.00150 

.00201 

.00120 

.00268 

.OO2U8 

-.00120 

-.00158 

.00316 

.0052 


The lead thermocouple was located 0,25^ centimeter from the surface 
and the convergence criterion was 10~^, 












TABLE IV.- COMPARISON OF CALCULATED HEAT RATES FOR 


CARBON-CARBON AND LI-900 FOR TRAJECTORY D^ 


Time, 

sec 

Known heat 
2 

rate, W/m 

Cedculated heat 
2 

rate, W/m 

Error, 

W/m^ 

Error , 
percent 

100 

18 158. U 

18 158.31U9 

0.08512 

O.OOOl+T 

200 

38 586.6 

38 586.628U 

-. 0283 T 

-.OOOOT 

300 

90 T 92.0 

90 T 1 *T .1136 

I+I+.89 

.0I+9I+I+ 

Uoo 

226 980.0 

226 98O.T60U 

.T60I+ 

-.00031+ 

500 

3 Ut 2 T 9 .** 

3I+T 28O.OU69 

-. 61+69 

-.00186 

6oo 

3 T 9 056.6 

3 T 9 01+5.8638 

10 .T 1 + 

.00283 

TOO 

31 k 5IT.0 

3 Tl+ 501.30I+3 

15. TO 

.001+19 

800 

3 T 2 2 kl .2 

3 T 2 26T.3I+I+5 

-20.11+ 

-.0051+1 

900 

3 Tl+ 5IT.O 

3 T 1 + 558.3558 

- 1 + 1.36 

-.01101+ 

1000 

3 T 9 056.6 

3 T 9 068.6OT2 

-12.01 

-. 0031 T 

1080 

381 326.1+ 

381 331 . 8 T 02 

- 5 . 1 +TO 

-. 0011+3 

1100 

3 T 6 T86.8 

3 T 6 T8T. 21+26 

-.Ul+26 

-.00012 

1200 

283 T 25.0 

283 T26.TO2I+ 

- 1 .T 02 

-.OOOTO 

1300 

192 933.0 

192 933.2610 

-.2610 

-.00011+ 

lUoo 

12k 839.0 

12I+ 936.5560 

-. 9 T 56 

-.OT815 


The thermocouple vas located on the hackwall of the carbon, 
and the convergence criterion was lO”^. 


IT 






Figure 2.- Heating rat , as a function of time for 
trajectories A, B, and C. 
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1000 I- 



0 .2 .4 .6 .8 

Thermocouple depth, cm 


Fig\ire 3.- Effect of average error due to thermocouple depth 
for a convergence criterion of 10”^. 


Pressure, N/m (atm) 


101 325 



Time, sec 


Figure 1*.- Pressure as a function of time for trajectory B. 
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Convective heating rate, MW/m 



APPENDIX A 


i 

I 

' CALCULATING SURFACE HEAT FLUX FROM KNOWN SURFACE TIMPERATURES 


To predict the heating rate at a siirface node, the houndary condition 


^3T 

3x 


x=0 


= ■^let 


(2a) 


is used. Using an implicit formulation, the finite difference equations can 
be rearranged into a tridiagonaJ. system of the form 


a T’ , +bT'+cT' _+d=0 
m m+l mm m m-1 m 


( 7 ) 


To visualize the relationship of the coefficients used to solve for 


%et 


gether with the other unknown temperatures, the matrix formulation is examined. 


b^a^OOOO OO 0 0 


\et 



0 b_ a« 0 0 0 0 


To 


do 

2 2 


2 


2 

0 c* b„ a_ 0 0 0 


Ti 


d^ 

3 3 3 


3 


3 

o 

o 

o 

o' 

XT- 

(0 

tr 

0 

1 
1 
1 
1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o 


Ti 

+ 

^1. 



T* , 


• 

• 

• 

d , 

''n-1 "n-1 n-1 


n-1 


n-1 

0 0 0 c b 


T» 


d 

n n 


n 


n 


(Al) 

The Gauss elimination method, discussed in reference 15, is applied to 
solve the system of equations. The solution of this matrix gives the tempera- 
ture of each node in the system for the next future time step. 


t. /• 




2h 



APPENDIX B 


THEFMOPHYSICAL PROPERTIES OF MATERIALS 


This appendix contains supportive data (tables B-I to B-IIl) concerning 
the thermophysical properties of reusable surface insulation (RSI) and reusable 
carbon-carbon (RCC) materials used in this investigation. 
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TABLE B-I.- THERMOPHYSICAL PROPERTIES OF RSI 



















TABLE B-II.- THERMOPHYSICAL PROPERTIES OP RSI COATING 


[Density l666 kg/m^] 


Tonperatvire , K 

Thermal conductivity, 
W/(m*K) 

Specific heat , 
J/(kg‘K) 

117 

0.735 

628 

172 

.788 

711 

256 

. 81+3 

795 

39k 

.951 

900 

533 

I.0U5 

lOOU 

672 

1.131 

1088 

811 

1.218 

1192 

950 

1.297 

1355 

1089 

1.376 

1318 

1200 

I.U3I+ 

1360 

1228 

I.UU9 

1381 

1339 

1.506 

1U23 

1367 

1.528 

Ikkk 

1U22 

I.5U9 

1I+6U 

II+50 

I.56U 

1U77 

1533 

I.61I1 

1506 
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